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ABSTRACT 

We investigate the properties of r-mode oscillations of a slowly rotating neutron star with 
a solid crust, by taking account of the effects of the Coriolis force. For the modal analysis we 
employ three-component neutron star models that are composed of a fluid core, a solid crust 
and a surface fluid ocean. For the three-component models, we find that there exist two kinds of 
r-modes, that is, those confined in the surface fluid ocean and those confined in the fluid core, 
which are most important for the r-mode instability. The r-modes do not have any appreciable 
amplitudes in the solid crust if rotation rate of the star is sufficiently small. We find that the 
core r-modes are strongly affected by mode coupling with the crustal torsional (toroidal) modes 
and lose their simple properties of the eigenfunction and eigenfrequency as functions of the 
angular rotation velocity fl. This indicates that the extrapolation formula, which is obtained in 
the limit of — > 0, cannot be used to examine the r-mode instability of rapidly rotating neutron 
stars with a solid crust unless the effects of mode coupling with the crustal torsional modes are 
correctly taken into account. 

Subject headings: instabilities — stars: neutron — stars: oscillations — stars: rotation 

1. Introduction 

Since the discovery of the gravitational radiation driven instability of the r-modes by Andersson (1998) 
and Friedman & Morsink (1998), a large number of studies on the properties of r-modes and inertial modes 
of rotating stars have been done to prove their possible importance in astrophysics (for a recent review 
see, e.g., Friedman & Lockitch 1999). Although early investigations of the r-mode instability were mainly 
applied to young hot neutron stars (e.g., Lindblom, Owen, & Morsink 1998; Owen et al. 1998), recently 
some authors have also discussed possible roles of the r-mode instability in old and cool neutron stars with 
a solid crust and a magnetic field, as those found in low mass X-ray binaries (LMXBs) (e.g., Andersson, 
Kokkotas, & Stergioulas 1999; Bildsten & Ushomirsky 2000; Rezzolla, Lamb, & Shapiro 2000). For example, 
Bildsten & Ushomirsky (2000) have suggested that the r-mode instability can be largely weakened by the 
effects of viscous damping in the boundary layer at the interface between the solid crust and the fluid core. 
Rezzolla et al. (2000) have also suggested that amplitude of r-mode oscillations tends to reduce due to the 
coupling to the magnetic field of a star. 

It is well known that neutron stars with a solid crust show a rich spectrum of nonradial oscillation 
modes, for example, nonradial modes associated with the solid crust and those associated with the fluid 
core and the surface fluid ocean (McDermott et al. 1988). If the effects of rotation are taken into account. 
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it is quite common to find avoided crossings, as functions of f2, between various oscillation modes of the 
neutron star model with a solid crust (Lee & Strohmayer 1996). This is also the case for r-modes. When 
the angular rotation frequency f2 of the star is small, the r-mode frequency observed in the corotating frame 
may be given by 

"F(FTT)' 

where I' and m are the indices of the spherical harmonic function representing the dominant toroidal 
component of the displacement vector. On the other hand, the crustal toroidal mode frequency may be 
given by (e.g., Strohmayer 1991) 

cu.in) ^ iUtiO) + j^^^y (2) 

where ijJt(O) is the oscillation frequency of the mode at = 0. The frequencies of the two modes get close 
to each other at 

I'd' + 1) 

n„oss- ^ ' u;t{0). (3) 
m 

As shown by McDermott et al. (1988) and Lee & Strohmayer (1996), the oscillation frequency of the 
crustal toroidal mode that has no radial nodes of the eigenfunction is of order of uJt{0)/ ^/GM/ ~ 10~^ 
for = 2, and hence we have Across/ \J GMj ~ 10~^ x (/' + 1) for the modes with I' = m, where M 
and R are the mass and radius of the neutron star, and G is the gravitational constant. This suggests 
that the mode crossing can happen between the two modes at sufficiently small rotation rates. Almost all 
previous discussions on the r-mode instability, however, did not take account of the effects of mode crossing 
between the two different kinds of modes .0 In fact, the discussions have been based on the assumptions 
that the modal property of the r-modes is almost independent of the rotation frequency il and hence that 
an extrapolation formula, obtained in the limit of $7 ^ 0, is applicable to rapidly rotating neutron stars 
with 57/ a/ GAI I i?3 ~ 1. Obviously these two assumptions are not necessarily correct for the r-modes that 
experience mode coupling with other kinds of oscillation modes. 

In this paper, we examine the r-mode instability of neutron stars that have a solid crust, by considering 
the effects of the mode crossing between the r-modes and the crustal toroidal modes. The plan of this paper 
is as follows. In §2, we show the modal properties of the r-modes of neutron stars with a solid crust. In §3, 
we discuss the dissipation timescales of the r-modes. §4 is for discussions and conclusions. 



2. r-Modes of Neutron Stars with a Solid Crust 

Neutron star models we use in this paper are the same as those used in the modal analysis by 
McDermott et al. (1988). The models are taken from the evolutionary sequences for cooling neutron 
stars calculated by Richardson et al. (1982), where the envelope structure is constructed by following 
Gudmundsson, Pethick & Epstein (1983). These models are composed of a fluid core, a solid crust and a 
surface fluid ocean, and the interior temperature is finite and is not constant as a function of the radial 
distance r. The models are not isentropic and the Schwarzschild discriminant A has finite values in the 
interior (see, e.g., Yoshida & Lee 2000b). 



^After our submission of this paper, a paper by Levin & Ushomirsky (2000) concerning a similar problem appeared in a 
preprint server "astro-ph" . They treated r-modes in a constant density neutron star with a constant shear modulus crust. 
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The method of caleulation of nonradial osciUations of rotating neutron stars with a soUd crust is the 
same as that used by Lee & Strohmayer (1996), in which the Cowhng approximation is employed. The 
terms due to the Coriohs force are included in the perturbation equations but no effects of the centrifugal 
force are considered. This approximation may be justified for low frequency modes satisfying |20/a;| > 1 
and Q"^ / {GAl / R^) ^ 1, where lo is the oscillation frequency observed in the corotating frame of the star 
(Unno et al. 1989). The eigenfunctions are expanded in terms of spherical harmonic functions Y[^{9,Lp) 
with different values of I for a given m. Here spherical polar coordinates (r, 9, ip) are used. For example, the 
Lagrangian displacement vector, is expanded as 

oo 

ir = rY, 5;(r)rr(0,^)e-*, (4) 

l>\m\ 

where I = \m\ + 2k and /' = i + 1 for even modes and I = \m\+2k + 1 and T = ^ — 1 for odd modes 
where fc = 0, 1, 2 • • • (Lee & Strohmayer 1996, see also Lee & Saio 1986). The symbol a denotes the 
oscillation frequency observed in an inertial frame, and we have lo — a + mVl. In this paper, to obtain a 
good convergence of the eigenfunctions and cigcnfrequencies at a given fJ, we keep an enough number of 
terms in the series expansion of the eigenfunctions. 

We computed frequency spectra of r-modes, inertial modes, and crustal toroidal modes for the 
three-component neutron star models called NS05T7, NS05T8, and NS13T8 (see, McDermott et al. 1988). 
The obtained mode spectra for the three models are qualitatively the same. In this paper, we therefore 
show the mode spectrum for the model NS13T8, because this is the most compact model among the three. 
The mass M and the radius R of the model NS13T8 are M = 1.326 Mq and R = 7.853 km, respectively, 
and the central temperature is Tc = 1.05 x 10^ K. 

In Figure 1, scaled eigenfrequencies k = w/il of r-modes, inertial modes, and a crustal toroidal mode of 
the model NS13T8 are given as functions of = n/^/GM/R^ for m = 2, where u) is the frequency observed 
in the corotating frame of the star. Here only the fundamental r-modes with V = m = 2 are considered 
since they are most important for the r-mode instability of neutron stars (see Lockitch & Friedman 1999; 
Yoshida & Lee 2000a, 2000b). The crustal toroidal mode shown in Figure 1 is the fundamental mode with 
no radial nodes of the eigenfunction and has the smallest oscillation frequency of the mode of this kind for 
a given V ai VI = Q (sec McDermott et al. 1988). Note that at O it is practically impossible to correctly 
calculate rotational modes because of their coupling with high overtone (?-modes. In this figure, we have 
used the notation given by iir^ and ^/r* for the r-modes, idn for the crustal toroidal modes, and i^i for the 
inertial modes, where V denotes the index of the spherical harmonic fimction associated with the dominant 
toroidal component of the displacement vector, n is the number of radial nodes of the eigenfunction, and Iq 
is the number introduced to classify the inertial modes (see, e.g., Yoshida & Lee 2000a). Note that the mode 
classification has been done at sufficiently small values of O. Since there appear surface r-modes and core 
r-modes for the three-component neutron star models, we have also introduced the superscripts s and c to 
distinguish the two kinds of r-modes. The oscillation energy of the core r-mode is predominantly confined 
in the fluid core, and that of the surface r-mode is confined in the surface ocean. A typical eigenfunction 
of a surface r-mode is given in Lee & Strohmayer (1996). The modes in the flgure are all odd modes (with 
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l' = \m\ for the toroidal modes) and they are retrograde modes propagating in the opposite direction to 
that of steUar rotation. 

Figure 1 clearly shows that the mode crossing between the 2''o"™^ode and the 2io-mode, as predicted 
by equations (0) to (^), leads to an avoided crossing of their frequencies as functions of At the avoided 
crossing, the modal properties of the two modes are exchanged. This is understood by examining the 
eigenfunctions of the two modes before and after the crossing, which are indicated by the filled circles in 
Figure 1, where the labels (a), (c), and (d) attaching to the filled circles are corresponding to the panels 
(a), (c), and (d) in Figure 2. Note that the 2io-niode at fj = 4 x 10"'^, which should carry the label (b), is 
not shown in Figure 1, because k ^ 1. In Figure 2, the expansion coefficients S3, H3, and iT2 of the two 
modes with Z' = m = 2atj7 = 4x 10^'^ in panels (a) and (b) and at = 4 x 10^^ in panels (c) and (d) 
are given as functions of the fractional radius r/R, where the eigenfunctions are normalized as iTii{r{,c) = 1 
with Tfjc being the bottom of the solid crust of the model. When Cl is sufficiently small, the 2'rQ-Tnode 
in panel (a) and the 2tQ-'niode in panel (b) have their typical eigenfunctions and the influence of other 
modes on the eigenfunction is negligible. Through the avoided crossing, however, the modal properties 
of the two modes are exchanged along the frequency curves as Cl increases. As shown by Figure 2, the 
mode in panel (d), which is corresponding to the mode (d) on the frequency curve of 2^0 in Figure 1, has 
in the fluid core the eigenfunction characteristic of the 2?'o"™'3de. On the other hand, the mode in panel 
(c), which corresponds to the mode (c) on the frequency curve of 2^0^ the crust the eigenfunction 

characteristic of the 2io-niode. However, this avoided crossing between the two modes is not sharp in the 
sense that the eigenfunction of each of the two modes suffers the contamination of the other mode after the 
avoided crossing. For example, the mode (d), which may be regarded as a core r-mode, has also oscillation 
amplitudes in the solid crust. This is in a sharp contrast with avoided crossings, for example, between 2?'o 
and 2^0 between 51 and 2''o' which the mode properties are almost completely exchanged through 
the crossings. Because of the avoided crossing between the 2^0"™^*-"^^ ^^"^ crustal toroidal modes, the 
2rg-mode will lose its simple modal property of the eigenfrequency and the eigenfunction as a function of 
Cl. We note that the fundamental r-mode with Z' = m of fluid stars has such a simple modal property as a 
function of n (e.g., Yoshida & Lee 2000b; Karino et al. 2000). 

Avoided crossing between the j/rg and ido modes with /' m may be common for neutron stars with 
a solid crust. To show this, we calculate the j/rg modes with I' = m = 2 for n ^ 1 polytropic neutron star 
models in which a solid crust is artificially embedded and isentropic structure is assumed. We assume the 
same physical parameters, such as the mass, radius, and the crust thickness, as those of the model NS13T8. 
For the shear modulus fiQ, we assume that the ratio fj,o/p^^^ is constant in the crust where p is the mass 
density (see e.g., McDermott et al. 1988). In Table 1, we tabulated the values of yL^j p^l"^ in the crust and 
the frequencies ijJt(G) of the fundamental 2^0 mode at = 0. We note that Wt(0) of the polytropic model is 
approximately proportional to the value of \/ /io/p^/^- Figure 3 gives the scaled frequencies k = w/fi of the 
2rg-mode as functions of fi for the polytropic models. The bends of the frequency curves of the r-modes as 

increases are caused by the avoided crossing with the crustal toroidal modes 2^07 as found in Fig.l for 
a realistic neutron star model. As Figure 3 indicates the avoided crossing between the 2T%-"ci\ode and the 
2io-mode is a common phenomenon and occurs at around ^cross^ which is determined by wt(0) for given /' 
and m (see equation (^)). The apparent difference in behavior of the r-mode frequencies as a function of il 
between the three polytropic models is therefore attributable to the difference in the location Across of the 
avoided crossing, that is, to the difference in the property pi^j p^^'^ of the solid crust. 

Here, it is worth while to note that, as shown by Yoshida & Lee (2000b), the deviation of the model 
from isentropic structure is not important for the fundamental r-modes with V = m and inertial modes. 
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if the Schwarzschild discriminant A of the nonisentropic model has sufficiently small absolute values as is 
expected in neutron stars. As a matter of fact, we can obtain similar results to those shown in Figure 1, 
even if we artificially set A = in the model NS13T8. 



3. Dissipation Timescales of Core j-Modes 

For comparison with other studies on the r-mode instability, let us derive an extrapolation formula 
for the damping timescale of the fundamental core r-mode with /' = m = 2 by calculating the r-modes at 
sufficiently small values of Q, for which the mode coupling with the crustal toroidal mode is negligibly weak. 
Then, using the same shear and bulk viscosity coefficients as those employed in Yoshida & Lee (2000a) 
and taking account of the viscous boundary layer (VBL) damping effects following the formulation given in 
Bildsten & Ushomirsky (2000) (see also Andersson et al (2000) for a correction), we may approximately 
express the formula as (see e.g., Lindblom et al. 1998) 

1 1 /io« 1 f Y f 



1 /lOS K\ / 02 \'^^ If 



TvBL \ Tc J \nGpJ fj^2 \ttGp 



(7) 



where p is the average density of the star, and Tc is the central temperature. The same viscosity coefficients 
are used both in the fluid and solid regions, and no effects of superfluidity in the core and in the inner 
crust are considered. Here, the first, second, third and fourth terms in the right-hand side of equation (Q) 
are contributions from the shear viscosity, the bulk viscosity, the shear viscosity in VBL and the current 
quadrupole radiation, respectively. The factors containing the ratio Tc/10® K have been introduced in 
equation (^ to extrapolate the damping rates with respect to the central temperature Tc- 

In Table 2, the various dissipation timescales f of the 2'"o modes with V ^ in ~ 2 are tabulated for the 
model NS13T8 and a simple n = 1 polytropic model with a static solid crust. For the polytropic model, 
we employ M = 1.4Mq and R — 12.5km and assume that a solid crust forms in the density region of 
p < 1.5 X 10^"* g/crn^ . Note that no wave propagation is assumed in the static crust of the polytropic model 
employed in this section. Bildsten & Ushomirsky (2000) (see also Anderson et al. 2000, and Rieutord 2000 
for more detailed treatment of VBL) used the similar polytropic model with a static solid crust to estimate 
the VBL damping effects on the r-modes. To compute the r-modes of the simple polytropic model, the 
method given in Yoshida & Lee (2000a), in which effects of the centrifugal force and perturbations of the 
gravitational potential are taken into account, has been employed. As shown by Table 2, all the dissipative 
timescales, except that due to the bulk viscosity, have similar values both for the simple polytropic model 
and for the model NS13T8. We note that although the destabilization due to the current quadrupole 
radiation for the model NS13T8 is by one order of magnitude stronger than that for the polytropic model, 
the VBL damping timescales, which are most important among the damping mechanisms of the r-modes, 
are almost the same for the two models. We consider that the difference in tj_2 between the two models is 
caused by the difference in the compactness GM / [c^R) of the models, where c denotes velocity of light. 
The huge difference in the damping timescales associated with the bulk viscosity between the two models 
may come from the difference in the treatment of the modes in the crust. Although we count the damping 
contributions in the crust for the model NS13T8, we do not count them for the polytropic model because 
we have assumed that the crust is static and allows no oscillations. In addition to this, the model NS13T8 
in an evolutionary sequence of cooling neutron stars has the crustal temperature that is higher by a factor 
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7 than the core temperature. Since the bulk viscosity coefficient is proportional to T^, the contribution of 
the bulk viscosity in the crust is largely enhanced for the damping of the /'-modes. 



4. Discussion and Conclusion 

It is well known that retrograde oscillations whose frequency uj satisfies the condition < uj/fl < m are 
unstable to the gravitational radiation reaction (Friedman & Schutz 1978). Using equation (||), we can show 
that the retrograde crustal toroidal modes satisfy this condition and become unstable to the gravitational 
radiation reaction when 

where ujt{0) denotes the oscillation frequency of the mode at 57 = 0. This instability can happen at slow 
rotation rates of the star if ujt{0) is sufficiently small. Note that all crustal toroidal modes shown in Figure 
1 are unstable to the gravitational radiation reaction. Our preliminary stability calculation, however, shows 
that the instability of the crustal toroidal mode is very weak and may not be important. 

In this paper, we have investigated the properties of r- modes of slowly rotating neutron stars with a 
solid crust. We have found that the mode property of the core tq mode with I' = m = 2, which is most 
important for the r-mode instability of neutron stars, is strongly affected by mode coupling with the crustal 
2to modes as D, increases. This means that we cannot assume a simple mode property of the core tq modes 
with I' — m a,s a. function of fl. To discuss the r-mode instability of rotating neutron stars it has been 
common to derive an extrapolation formula for the damping timescale of the r-modes (e.g., Lindblom et al. 
1998). The extrapolation is usually carried out about the interior temperature T and the rotation frequency 
il. The extrapolation with respect to the temperature may be justified if the stars have an isothermal 
structure characterized, for example, by the central temperature Tc- On the other hand, since the mode 
properties of the r-mode are well known only in the limit of £7 — ^ 0, the damping timescales are calculated 
in this limit and are extrapolated to rapidly rotating neutron stars with 0^1. This process of the 
extrapolation with respect to is justified if the mode frequency w is well approximated by 2mn/l'{l' + 1) 
and the eigenfunction has the dominant toroidal component, independent of O. We know that the above 
two assumptions are satisfied for the fundamental r-modes with V — m for fluid stars (Yoshida & Lee 2000b; 
Yoshida et al. 2000; Karino et al. 2000). In the case of neutron stars with a solid crust, however, the two 
assumptions are not satisfied except when f2 ~ and no mode couplings are expected. For example, when 
the r-mode is near the avoided crossing with the crustal 2^0 niode, the eigenfunctions and eigenfrequencies 
of the r-mode behave peculiarly and change rapidly as functions of 17. This makes it difficult to predict 
reliably the evolution of a neutron star driven by the r-mode instability through the avoided crossing by 
simply extrapolating the damping timescales obtained at small f7. In addition to this, the eigenfunctions 
of the r-modes are contaminated by the crustal torsional modes at large 17, and the extrapolated damping 
timescales are not necessarily accurate enough to give the useful estimations. These difficulties mean that 
the extrapolation formula of the damping timescales of the r-modes, derived in the limit of 17 ^ 0, is not 
applicable to the discussion on the r-mode instability of the rapidly rotating neutron stars with a solid crust 
unless the effects of the avoided crossings with and the contamination of the eigenfunctions by the crustal 
torsional modes are correctly taken into account. In this sense, discussions based on the extrapolation 
formula of the damping timescales are not fully justified for neutron stars with a solid crust. 
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Table 1. Shear modulus /xq and eigenfrequency u)t{0) of the 2*0 
mode for the polytropic model with a crust at O = 0. 



Model 




t.;t(0)/(GM/i?3)i/2 


model 1 


1.0 X 10" 


9.9 X 10- 


-3 


model 2 


1.0 X 10^^ 


3.1 X 10" 


-2 


model 3 


1.0 X 10^^ 


9.9 X 10" 


-2 



Table 2. Dissipative timescales f of core ro modes with /' = m = 2 for neutron star models. 



Model 




rs{s) 


fvBL{s) 




NS13T8 


1.3 X 10'' 


1.6 X IQB 


3.7 X 10 


-2.2 X IQ-i 


Polytrope with a static crust 


1.7 X 10^^ 


1.9 X 10« 


8.9 X 10 


-4.4 X 10° 



1 



.8 



.6 



.4 



.2 

.01 .02 .03 .04 .05 

Q/(G M/r')'/' 

Fig. 1. — Scaled frequencies k = ui/Sl of the r-modes, inertial modes and the crustal torsional (toroidal) 
mode of the model NS13T8 arc plotted as functions of = ^/ [GM / E?)^/^ for m = 2, wlioro the frequency 
curves are attached by the labels indicating their classification at small values of f7 (see the text for the 
detail of the classification). The modes depicted are retrograde odd modes. The r-modes and the crustal 
toroidal mode have the dominant toroidal component associated with Z' = m = 2 at small values of f2. Note 
that the eigenvalue of the r-modes is given by k = 2m/{l'{l' + 1)) in the limit of f2 — »■ 0. The filled circles 
indicate the modes whose eigenfunctions are shown in Figure 2. 
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0.2 0.4 0.6 0.8 10 0.2 0.4 0.6 0.8 1 

r/R r/R 

Fig. 2. — Expansion coefficients ^3 (solid curve), (dashed curve), and iT^ (dotted curve) of the core 
r-modes and the crustal toroidal modes with T = m = 2 for the model NS13T8 are given as functions of r/R 
at n = 4 X 10~^ in panels (a) and (b), and at = 4 x 10~^ in panels (c) and (d), where the normalization of 
the eigenfunctions are given as iT2{rbc) = 1 with r(,c being the bottom of the solid crust. The eigenfunctions 
shown in panels (a) and (b) are those of the 2fQ-Tiiode and the 2io-iiiode, respectively. The locations of the 
modes are shown in the k — ^ plane by the filled circles in Figure 1. 
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Fig. 3. — Scaled frequencies k = u/fl for the fundamental 2?'o modes with I' = m = 2 are plotted as functions 
of = ^/{GM/R^y^^ for the simple polytropic model with a sohd crust. The three models have different 
values of /xq/p*^^ that is assumed to be constant in the crust (see, Table 1). 



